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摘 要 传统 的 有 中 介 的 调节 (mediated moderation, meMO ) 模型 关于 误差 方差 齐 性 的 假设 
经 常 被 违背 ， 应 用 研究 中 也 缺乏 测量 meMO 效应 大 小 的 指标 。 对 于 单 层 数据 ， 本 文 借助 于 


两 层 建 模 的 思想 , 提出 了 一 种 可 用 于 处 理 方差 


非 齐 性 的 两 层 有 中 介 的 调节 (2meMO ) 模型 ; 


给 出 了 用 于 测量 meMO 分 析 中 总 调节 效应 、 直 接 调节 效应 和 有 中 介 调 节 效 应 大 小 的 效应 量 


通过 Monte Carlo 模拟 研究 ， 比 较 了 meMO 玫 


I 2meMO 模型 在 参数 和 效应 量 估计 上 的 表 


并 通过 实际 案例 解释 了 2meMO 模型 的 应 用 以 及 效应 量 的 计算 和 解释 。 
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现 。 


在 很 多 领域 的 实际 研究 中 ， 中 介 (mediation ) 效应 分 析 和 调节 (moderation ) 效应 分 析 


poll 


是 探讨 自 变量 X MARA 
统 的 理论 框架 下 , 将 中 介 效 应 与 调节 效应 结合 


TH 


了 之 间 复 杂 关 系 常用 的 两 类 方法 。 随 着 研究 问题 的 深入 ， 在 更 系 


， 探 讨 变量 之 间 更 为 复杂 关系 的 研究 方法 越 来 


越 受 到 重视 ( Kwan & Chan, 2018; 叶 宝 娟 , WERE, 2013)。 其 中 ， 有 中 介 的 调节 (mediated 


moderation， 简 称 meMO) 模型 在 探讨 调节 变量 作用 机 制 方面 ， 提 供 了 强 有 力 的 分 析 方 法 。 


meMO 模型 分 析 的 重点 是 探讨 调节 变量 Z 对 自 变 量 (X) 与 因 变 量 (7) 之 间 关 系 的 作 


机 制 ， 如 果 Z 对 并 和 了 之 间 关 系 的 调节 通过 


起 作用 ， 则 称 2 通过 M 对 和 和 了 之 间 关 系 间 接 的 调节 为 有 中 介 的 调节 (meMO )。 


中 介 的 调节 变量 (M) (mediating moderator) 
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理论 的 假设 , meMO 模型 近年 来 在 心理 学 研究 


基于 不 同 


得 到 了 广泛 的 应 用 (如 , Baron & Kenny, 1986; 


Liu et al., 2012; Muller et al., 2005; 王 玲 晓 等 , 2019; 杨 文 圣 等 , 2019; 杨 逸 群 等 , 2020)。 在 方 


法 研究 领域 ， 也 有 许多 古 


究 者 就 meMO 模型 的 分 析 步 又 和 检验 方法 进行 了 较为 详尽 的 介绍 


和 总 结 (Hayes, 2018; Kwan, & Chan, 2018; Ng et al., 2019; 叶 宝 娟 , WERE, 2013)。 然 而 ， 在 


实际 应 用 中 meMO 模型 仍 至 少 面临 以 下 三 个 方面 的 问题 。 

首先 ， 对 于 meMO 模型 ， 由 
， 或 者 研究 中 所 选 有 
调 统计 显著 性 的 问题 , 越 来 越 多 的 期 


其 效应 


E 
tee Pit ap 


的 替代 性 


里 的 meMO 效应 量 的 指标 ， 应 用 研究 中 # 
指标 与 meMO 研究 感 兴趣 的 问题 不 匹配 。 针 对 研 
站 要 求 研究 者 提供 效应 量 和 
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区 间 估 计 的 结 


=j 


Rs 
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RW, American Education Research Association, 2006; Cumming, 2014; Funder et al., 2013; 


Rozeboom, 1960; Wilkinson & American Psychological Association Task Force on Statistical 


Inference, 1999), 


国内 的 《心理 学 报 》 和 《心理 也 有 同 档 


的 要 求 。 然 而 ， 到 目 


前 为 止 ， 研 究 者 没有 结合 meMO 问题 的 特点 对 其 效应 量 进行 针对 性 的 测量 和 解释 。 有 研究 


者 借鉴 其 他 类 似 的 效应 量 作为 替代 品 ， 如 采用 
调节 路 径 后 ， 因 
应 量 定义 的 本 质 来 看 ， 效 应 量 是 反应 某 种 现象 (如 
应 量 应 直接 针对 所 感 兴趣 的 研究 问题 (Kelley & Preacher, 2012) 。 因 此 ， 吸 待 解决 的 问 
如 何 围绕 meMO 分 析 感 兴趣 的 “Z 通过 M 调节 和 和 了 之 | 
HR 
其 次 ， 有 中 介 的 调节 (meMO) 和 有 i 


EE 点 不 同 ， 但 是 在 实际 应 月 


然 解释 的 侧 习 


量 指标 ， 对 meMO 4474 


Lambert, 2007; Liu et al., 2012; Muller et al., 2005). Hayes (2018) 


变量 变异 解释 


| 


+ mt 


的 有 中 介 的 调节 效应 大 小 ; 


常用 的 AR? 描述 增加 了 有 中 介 的 
率 的 增加 (Aiken etal., 1991)， 但 这 些 指标 并 不 合适 。 因 为 从 效 


啊 程 度 ) 大 小 的 数量 化 表示 ,合适 的 效 


日 


LIE 


司 的 关系 ”问题 ， 构 建 合理 的 效应 


关节 的 中 介 (moderated mediation， 简 称 moME) 虽 
MIA. Baron 和 Kenny(1986) 最 早 提 
出 了 meMO 和 moME 的 概念 ， 很 多 研究 随后 就 这 两 类 分 析 进 行 了 详细 的 阐述 (Edwards & 
概念 模型 和 统计 模型 
对 meMO 和 moME 进行 了 详细 的 分 类 ， 并 提供 了 PROCESS 分 析 软 件 。 然 而 ， 对 于 应 ) 
究 者 来 讲 清晰 区 分 两 类 模型 并 不 容易 ， 丰 
moME (EÈ meMO) 的 叙述 ， 而 在 统计 模型 结果 的 解释 上 却 采 上 月 


$ 


Jø 
会 出 现 理论 假设 中 变量 之 间 的 关系 是 


月 了 meMO (或 moME) 的 解 


释 ( 如 , Gonzalez-Mulé et al., 2016; Liang & Chi, 2013; Park & Searcy, 2012; Tang, 2016; Tang & 


Naumann, 2016; Zacher et al., 2012; 杨 文 圣 等 , 2019)。 究 其 


meMO 与 moME 对 应 的 统计 模型 相同 ， 导 致 丰 


$e 


pii 


型 等 价 ，meMO-II 模型 与 


,常用 的 meMO 模型 分 为 meMO-I 型 和 meMO-I 型 下 


AAA — 


a 


EAN TER E, 
究 者 对 meMO 和 moME FE SB. RRK 
类 , meMO-I 模型 与 第 一 阶段 moME 


阶段 meMO 模型 等 价 (Kwan & Chan, 2018)。 本 文 关 注 的 问 


I 


题 是 如 何 从 概念 模型 出 发 ， 在 统计 模型 的 表述 上 突出 meMO 模型 关注 的 变量 间 关 系 以 及 调 
节 变 量 的 作用 机 制 ， 定 义 合理 的 效应 量 ， 以 在 结果 解释 上 区 分 meMO 和 moME 两 类 模型 。 


最 后 ， meMO 模型 假设 误差 方差 齐 性 ,回归 系数 显著 性 的 检验 也 基于 这 一 假设 , 然而 实 


际 应 用 中 这 一 假设 常常 不 成 立 ， 尤 其 是 包含 调节 效应 的 模型 (Aguinis & Pierce, 1998; DeShon 


& Alexander, 1996)。 例 如 ， 对 于 含有 交互 作用 的 回归 模型 


Y; = by + bX; + bZ; + bsXiZi+ £i, (1) 
式 (1) Y XA Z 分 别 表 示 因 变量 、 自 变量 和 调节 变量 ; b 表示 非 标 准 化 的 回归 系数 ; 
误差 e:~N(0,0?)， i=1,2,. nmn 表示 样本 量 。 方 差 齐 性 的 假设 条 件 ， 意 味 着 对 于 自 变 量 
X 取 不 同 的 值 刀 时， 误差 方差 相等 。 然 而 ， 含 有 显著 交互 作用 的 模型 很 难 满 足 方差 齐 性 的 前 
ERR, 由 此 会 导致 第 二 类 错误 的 增 大 和 检验 力 的 下 降 ， 当 交互 作用 项 对 应 的 回归 系数 显著 
不 为 零 时 ， 其 参数 估计 对 应 的 标准 误 有 偏 ， 从 而 导致 推断 统计 的 结果 (包括 效应 量 ) 出 现 偏 


差 (Aguinis & Pierce, 1998; DeShon & Alexander, 1996) 。 本 研究 拟 对 传统 的 meMO 模型 进 
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S 行 拓展 ， 在 此 基础 上 定义 更 合理 的 测量 meMO 效应 量 的 指标 ， 同 时 也 解决 误差 依赖 于 自 变 
© 量 的 方差 非 齐 性 问题 。 

= 为 了 处 理 调节 效应 分 析 容 易 导致 的 方差 非 齐 性 的 问题 ，Yuan 等 人 (2014) 采 用 两 层 模型 
= 扩展 了 传统 的 调节 效应 回归 模型 ， 并 据 此 定义 了 测量 调节 效应 大 小 的 指标 。 虽 然 Yuan 等 人 
> (2014) 没 有 涉及 meMO 模型 ， 但 他 们 提出 的 两 层 建 模 的 思想 和 效应 量 定义 的 出 发 点 对 调节 


效应 拓展 到 更 复杂 的 模型 是 有 启发 的 。 本 研究 将 借鉴 两 层 模型 定义 的 思路 , 将 其 拓展 到 复杂 
的 有 中 介 的 调节 模型 的 情景 ， 并 基于 此 定义 既 适 用 于 方差 齐 性 条 件 ， 又 适用 于 方差 非 齐 性 条 
件 的 meMO 效应 量 指标 。 

综 上 ， 本 研究 主要 解决 以 下 问题 : 首先 从 有 中 介 的 调节 的 概念 模型 出 发 ， 将 Yuan 等 人 
(2014) 定 义 的 两 层 调节 回归 模型 拓展 到 两 层 有 中 介 的 调节 (2meMO) 模型 ; 其 次 通过 对 自 变 
me 了 与 因 变量 了 之 间 关 系 变异 来 源 的 分 解 ， 给 出 测量 meMO 效应 量 大 小 的 指标 ， 然 后 通过 
模拟 研究 评估 模型 参数 以 及 效应 量 估计 的 精度 ， 并 通过 与 传统 meMO 比较 ， 验 证 拓展 模型 
的 合理 性 ， 同 时 采用 实际 数据 案例 和 Mplus(Muthén & Muthén, 2017) 语 名 介绍 2meMO 模型 
的 应 用 以 及 效应 量 的 估计 和 解释 , 以 方便 应 用 者 使 用 ; 最 后 指出 方法 的 可 拓展 性 以 及 效应 量 
指标 的 优 缺 点 和 实际 应 用 中 应 注意 的 问题 。 


2 有 中 介 的 调节 模型 的 拓展 


虽然 meMO 模型 
角度 ，meMO-I 模型 更 有 意义 ， 有 很 多 下 
义 ， 应 用 中 不 应 被 关注 


I 可 以 分 为 meMO-I 型 和 meMO-IH 型 ,但 从 可 解释 性 和 变量 具体 意义 的 
FAUNA meMO-I 模型 中 对 应 的 变量 没有 实际 意 


E(Edwards & Lambert, 2007; Hayes & Preacher, 2013)。 因 此 ， 本 文 重点 


讨论 更 具 解 释 性 和 理论 意义 的 meMO-I 模型 (Kwan & Chan, 2018; Hayes, 2018)， 其 对 应 的 概 
图 1 中 如 果 Z 调 节 了 X 一 了 间 的 关系 ， 且 Z 在 X 了 Y 上 的 调节 效应 被 变量 M 
所 中 介 ， 这 时 就 产生 了 有 中 介 的 调节 《〈 粗 线条 标 出 的 间接 路 径 ) ， 变量 M 称 为 中 介 的 调节 
，M 可 能 只 是 部 分 中 介 了 Z 的 调节 效应 ， 所 以 Z EX 一 了 Y 上 的 直接 路 径 被 保 


念 模型 见 图 1 。 


变量 。 典 型 


。 图 1 的 概念 模型 只 是 各 变量 
能 够 反应 概念 模型 特征 的 统计 模型 ， 


之 间 关 系 的 理论 表达 ,不 能 直接 进行 分 析 。 只 有 将 其 转换 为 
才能 进行 进一步 的 数据 分 析 。 


图 1 有 中 介 的 调节 效应 的 概念 模型 


本 文 借助 于 Yuan 等 人 (2014) 的 思想 ， 采 用 两 层 建 模 的 思路 构建 图 1 对 应 的 统计 分 析 模 
型 。 图 1 的 概念 模型 重点 关注 针对 了 影响 程度 如 何 受到 调节 变量 的 影响 ， 即 对 于 不 同 的 个 
影响 程度 如 何 随 着 其 在 调节 变量 上 取 值 的 不 同 而 变化 。 第 一 层 模 型 描述 meMO 


ti, Xt YH 


分 析 中 自 变量 


响 ， 以 及 调节 变量 Z 和 中 介 的 调节 变量 M 之 
系 的 变异 ， 在 考虑 了 模型 中 
有 解释 的 变异 ， 


层 1: 


层 2: 


了 对 因 变 量 了 影响 程度 的 变化 ， 第 二 层 描述 Z 和 M 对 对 和 了 之 间 关 系 的 影 


因此 在 模型 ， 


间 的 关系 。 值 得 注意 的 是 ， 对 于 蕊 和 了 之 间 关 


涉及 到 的 调节 变量 Z 和 中 介 的 调节 变 量 M 外 ， 可 能 还 会 存在 没 


增加 X 一 了 影响 


的 随机 误差 项 。2meMO 模型 可 表示 为 : 


Y; = di + ciXi + £i (2) 
di = Yao + Ya Mi + YazZi + Uai (3) 
ci = Veo + YeiMi + YezZi + Uci (4) 
Mi = Ymo + YumiZi + emi (5) 
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方程 (2) - (4) H 
方程 (5) Hi 


量 了 和 中 介 


Z Fl M 
Ue, 服从 正 态 分 
var(uci) = 02。。emi， eyi 和 uct 是 相互 六 
但 是 实际 数据 并 不 是 嵌 套 的 两 层 数 据 , 在 参数 估 
另外 ， 在 模型 的 定义 中 ， 


问题 ( 
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为 随机 


XIZ 


响 (除了 对 外 直接 指向 因 


R T É 


BREE X, 


调节 变量 Z 和 中 


述 了 调节 变量 Z 和 中 


的 调 
均值 


三 | 
Fe 


节 变 量 M 的 
心 化 的 变量 ， 即 
布 ， 且 BE(eri) = E(emi) = E(u) = 0, var(eyi) = Oey » var(eyi) = Oem ? 


mite, 62 


口 


JAE 


介 的 调节 变量 M 与 因 
介 的 调节 变量 M 之 间 的 关系 。 


E; Yar» Yaz» Yoo Yer? Ycz 和 ymi 为 回 


RIA AH 


计 上 没 法 区 


MARNIE 


了 详细 讨论 这 


层 1 


问题 


变量 了 


A Yao ty mot 


型 虽然 用 两 


分 ei 和 和 wg; 9 因 


过 间 的 关系 ; 
别 表示 因 变 
HR% Bis xX, 


E(Z;) = E(X) = E(M;) =0。 残 差 项 eri = ei + Uqi» emi ™ 


RMT, 
此 将 其 合并 为 err。 


层 统 计 模 型 可 用 图 


许 c; 存 在 随机 变异 ， 以 此 来 处 理 误差 方差 随 X; 变 化 的 方差 非 齐 性 
Ye 2meMO 模型 对 应 的 两 


2 表示 。 


图 2 拓展 的 有 中 介 的 调节 〈2meMO) 分 析 对 应 的 分 层 统计 模型 


一 层 模 型 


系数 ; 第 二 层 模型 } 


和 MM 之 i 


司 的 关系 。 


具体 来 讲 , 方程 (3) 描述 了 模型 中 


LER 
E: 


变量 了 的 变相 


); 方程 (4) 描述 了 XY 对 了 的 


J 描述 meMO 重点 关注 和 解释 的 关系 是 针对 了 的 影响 ， 其 中 截 
步 解释 Z 和 M 对 随机 系数 di; 和 ci 的 解释 , 同时 第 二 层 模 型 还 定 
归 的 截 距 di 受到 的 影 
回归 系数 ci 受到 的 调节 


YX X (al 


Hd, FUR Kc; 


(指向 对 到 了 回归 线 上 的 变量 ), 其 中 yoo 表 示 了 对 对 了 的 主 效应 , ya HAR MX XM YER 


[E] 
Hy 


系 的 i 


节 ，ycz 表 示 Z 对 XY 和 了 之 间 关 系 的 直接 调节 ; 方程 (5) 


种 采 


两 层 定 义 模型 的 方式 突出 了 meMO 模型 中 解释 的 习 


如 何 被 调节 (直接 调节 或 间接 调节 ), 从 而 可 有 效 区 分 meMO 与 moME 所 关注 问题 的 差异 。 
通过 合并 模型 的 表达 进 


下 


述 了 


一 步 明 确 2meMO 模型 


5 


重点 是 和 了 之 间 的 关系 ( 


I 和 传统 meMO 模型 


M 对 Z 的 回 


归 。 这 


z 1) 


4 的 联系 和 区 别 。 


将 (3) 和 (4) 代入 到 D 中 ， 得 到 合并 模型 的 表达 式 为 : 


Y, = dj + GA + ei 


= Yao + Yar Mi + YazZi + Uqi + (Yeo + YeıMi + 


Vc2Zi + Uci)X; + £i 


= Yao + VarMi + YazZi + YeoXi + Ver XiMi + Ve2XiZi + UciXi + eyi, (6) 


再 将 等 式 〈5) 代入 到 等 式 (6 中 ， 得 到 : 


Y; = (Yao + Ya1Ymo) + (Yaz + YaiYm1)Zi + (Yeo + Yc1Ymo)Xi 


+(Yez + Ve1¥m1)XiZi + Var emi + VerXiemi + u 
= (Ya, + Ya, Ymo) + (Ya, + Ya, Ymı )Zi 


+[(yeo + Yc1Ymo) + (Vez + Ver Yu) Zi + (Yeiem 


ciXi + eri 


it Uci)]Xi + (Varemi + eyi). (7) 


需要 说 明 的 是 ， 相 比 传统 的 meMO 模 型 ， 方 程 (6) 和 《7) 多 了 uciXi 这 一 项 ， 此 项 可 以 


有 效 区 分 模型 解释 的 重点 是 有 ， 


介 的 调节 , 而 非 有 调节 的 中 介 (我 们 将 在 讨论 部 分 进一步 解 


释 meMO 和 moME 的 区 别 )。 另 外 ， 合 并 后 的 方程 (7) 中 隐 含 着 X 一 XM 和 XZ 一 XM 两 条 路 
径 以 及 对 两 条 路 径 系 数 的 限定 〈 见 图 2 和 图 3)， 即 限定 X 一 XM 的 路 径 系数 与 M 对 X 回 归 的 截 


如 下 : 从 方程 CD) 可 看 出 ，X 对 7 的 影响 为 yco + yeiywo， 包 括 Y 对 7 的 直接 影响 yco， 以 及 X 通 


过 XM 对 7 的 间接 影响 ycjymo， XM 对 7 的 影响 为 yci， 则 可 通过 限定 X 一 XM 为 ymo， 使 得 X = 


XM 一 Y 为 yeiymo; 同 理 ， XZ 对 7 的 影响 y。。 十 YciyYwi， 包括 XZ 对 7 的 直接 影响 yc>， 以 及 XZ 通 


过 XM 对 7 了 的 间接 影响 yexym1，XM 对 7 的 影响 为 ye1， 可 通过 限定 XZ 一 XM 为 ymi1 使 得 XZ > 


图 3 所 示 。 对 比 图 2 和 图 3， 图 2 的 解释 更 加 直接 。 


日 于 2meMO 分 析 所 对 应 的 合并 统计 模型 如 


的 残 差 方 差 可 以 随 X 的 变化 而 变化 (方差 非 齐 性 ), 传统 
型 允许 包含 随机 项 uci 不 仅 不 再 需要 误差 方差 齐 性 的 假设 ， 而 且 更 符合 实际 情况 。 因 为 除了 
模型 中 考虑 到 的 调节 变量 ,XR 和 7 之 间 的 关系 还 有 可 能 依赖 于 模型 中 没 包含 的 其 他 调节 变量 ， 
响 的 随机 误差 项 uci 也 合理 的 刻画 了 这 一 点 。 除 了 这 一 随机 路 径 ， 该 模型 与 传统 的 
meMO 对 应 的 统计 模型 (Grant & Berry, 2011; Liu et al., 2012) 相 比 ， 传 统 模型 中 XZ > XM 和 
X 一 XM 的 路 径 自 由 估计 , 而 在 我 们 推导 的 统计 模型 中 这 两 条 路 径 分 别 限定 为 与 2 对 M 回 归 的 
斜率 和 截 距 相 等 ,对 于 模型 的 自由 度 增加 了 2。 除了 随 书 


XK YS 


图 3 中 X 到 7 路 径 系 数 上 的 实心 圆 点 表示 其 对 应 的 回归 系数 含有 随机 效应 uc;， 用 以 描述 7 


meMO 模 型 中 不 包含 这 项 。 2meMO 模 


[效应 wei 外 ， 昌 然 2meMO 对 应 合并 统 


计 模 型 


的 其 他 部 分 与 Kwan 和 Chan(2018) 从 另 一 视角 提 


E 导 的 meMO 的 统计 模型 相同 ,但 两 者 


建构 和 解释 meMO 影 响 效应 的 思路 存在 本 质 差异 。2meMO 分 层 统计 模型 《图 2) 与 概念 模型 


的 对 应 关系 更 清晰 ， 模 型 中 关于 X 和 Y 之 间 关 系 变异 的 解释 也 更 易于 理解 ， 可 以 围绕 X 和 7 之 


间 关 系 变异 Co 的 解释 , 定义 调节 变量 Z 的 总 调节 效应 、 直 接 调 节 
而 且 这 一 定义 效应 量 的 思路 同样 适用 于 wci = 0 的 模型 。 


Xi emi 


效应 以 及 间接 调节 效应 ， 


图 3 拓展 的 有 中 介 的 调节 (2meMO) 分 析 对 应 的 合并 统计 模型 


3 meMO 效应 量 的 定义 


在 方程 (2) - (5〉 所 定义 的 2meMO 模 型 中 ,将 (5) 代入 (4) ， 可 以 将 X 对 了 总 的 影响 


效应 ci 表示 为 : 
Ci = Veo + Ye Mi + YezZi + Uci 


= Yeo + Yei(Ymo + YmiZi + emi) + YezZi + 


Uci 


= (Yeo + Ver ¥mo0) + (Vez + Ver ¥m1)Zi + Yeremi + Uci- (8) 
第 一 项 yco =F Vor Vmo EXAT YH Ziq, H Fx 中 yo 表示 对 7 的 直接 影 响 ， Yci1Ymo 表 不 XX 通 前 过 XM 对 
7 的 间接 影响 ， 这 一 项 与 调节 变量 Z 的 取 值 无 关 。 第 二 项 (yz + yciywi) 表 示 调 节 变 量 Z 对 X > 


Y 的 总 调节 效应 ， 其 中 yes 是 Z 对 X => 7Y 的 直接 调节 效应 ， Ver Yui eZ 


前 过 M 对 X 一 了 的 中 介 的 调 


节 效 应 (对 应 概念 模型 (图 2) 中 粗 体 线 表 示 的 Z 通 过 M 对 X 和 7 之 间 回 归 系 数 的 影响 ， 对 应 于 
合并 统计 模型 (图 3) 中 的 XZ 一 XM >Y) ，yerym1 即 为 mneMO 模 型 关注 的 有 中 介 的 调节 ， 
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记 为 ymemo。 推 断 统 计 中 将 中 介 调 节 效 应 的 检验 转换 为 检验 Ho:ymemo = 0 是 否 显著 的 问题 。 
实际 中 ， 经 常 通过 对 调节 变量 Z 取 给 定 的 值 z〈( 如 高 于 或 低 于 均值 1 个 标准 差 )， 计 算得 到 
Yc1Ym12， 来 解释 调节 变量 Z 取 不 同 值 时 ，2Z 通 过 中 介 变 量 M 对 X 和 7 之 间 关 系 的 影响 程度 〈 简 
单 meMO)。 类 似 于 调节 效应 中 简单 效应 差异 的 检验 ， 对 于 给 定 的 不 同调 节 效 应 值 Z1 = zi 和 
Z = Zz2， 可 得 到 中 介 调 节 效 应 差异 difiemo = YeryYmi(21 一 22)， 实际 中 也 可 对 Ho: difmemo = 
0 进行 检验 ， 以 回答 调节 变量 Z 取 不 同 值 z1 和 zs 时 ， 中 介 调 节 作 用 是 否 存 在 差异 。 

X 对 7 影响 程度 的 变异 ， 即 ci; 的 方差 为 : 


var(ci) = var (co 十 YeiMi + YezZi + Uci) 
= (yez + YeryYm1) var (Zi) + ya oe + ofe. (9) 
下 面 我 们 基于 ci 变异 的 不 同 来 源 来 定义 可 以 测量 meMO 大 小 的 效应 量 。 第 一 个 效应 量 


(Yez + Ye1¥m1)*var(Z;) 
(yez + Ve1¥m1)*var(Z;) + VA Oey + aĉe 
为 调节 变量 Z 对 回归 系数 cj; 总 变异 所 解释 的 比例 ， 可 以 用 来 测量 XY 对 7 影响 的 变异 中 ， 调 节 变 


量 Z 所 能 解释 的 总 变异 所 占 的 比例 。 效 应 量 


Yo tot = (10) 


Yevar(Z;) 
(yez + Ye1¥m1)2var(Zi) + VA Oey + Fic 
则 描述 了 变量 Z 的 直接 调节 效应 对 回归 系数 cj 总 变异 所 解释 的 比例 ， 可 以 用 来 测量 调节 变量 


Z 对 X 一 Y 路 径直 接 调 节 效 应 的 大 小 。 而 效应 量 


do air = 


(Ye1¥m1)?var(Z;) 
2 2: ed 2 (12) 
Ye + Yc1¥m1) var(Z;) 十 Yc1%em + Ove 


描述 了 变量 Z 通 过 中 介 的 调节 变量 M 的 间接 调节 效应 (也 叫做 有 中 介 的 调节 ) 对 回归 系数 cj 总 
变异 所 解释 的 比例 ， 可 以 用 来 衡量 调节 变量 Z 通 过 中 介 的 调节 变量 M 对 X 一 了 路 径 系 数 的 调 
节 效 应 的 大 小 。 值 得 注意 的 是 ， 基 于 变异 分 解 所 定义 的 效应 量 ， 由 于 直接 调节 效应 对 应 的 路 
径 系数 ycz 与 间接 调节 效应 对 应 的 路 径 系数 yeiywa 均 与 变量 Z 有 关 ， 两 者 之 间 不 具有 独立 性 。 
因此 , 通过 路 径 系数 定义 的 变异 解释 率 不 具有 可 加 性 , 其 直接 调节 效应 量 与 间接 调节 效应 量 


pmemo = do _ina = 


的 和 并 不 等 于 总 的 调节 效应 量 ， 即 : bao toc # modir 十 办 Mo inae 

当 方差 齐 性 假设 被 满足 时 , 回归 系数 ci 的 残 差 随机 部 分 为 零 , 即 ver = 0, 等 式 (2)- (5) 
中 的 模型 被 还 原 为 传统 的 有 中 介 的 调节 模型 (简称 为 mneMO), 上 面 定义 的 公式 (10) - (12) 
被 简化 为 以 下 公式 。 


总 调节 效应 ; 


D = — Wen * Yeas) varCD 
MOtot (Yoa + Ye1¥u1)2Var(Z;) + YE oy 


(13) 
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直接 的 调节 效应 : 


的 =- a 
Modir (Yez + yerymi) Var(Z2i) + ya ol 
中 介 的 调节 效应 : 
2 
N 4A = (Ye1¥u1)°Var(Z;) 
memo ~ MO_ind ` (15) 


(Yez 十 yeiywi)2var(Zi) + ya ofm 
这 里 ， 我 们 使 用 上 标 Im 来 区 分 定义 中 不 包含 的 X 对 了 未 解释 的 随机 残 差 。 


公式 (13) - C15) 同样 适用 于 2meMO 模型 ， 其 分 母 表 示 的 意义 是 在 不 考虑 X 一 了 影响 
系数 的 随机 误差 (ui;) 时 ， 调 节 变 量 Z 对 X 一 了 影响 系数 的 调节 。 其 包含 可 解释 的 调节 


(Yeo 十 YerYm1)*Var(Zi) 和 通过 中 介 的 调节 变量 M 间接 未 解释 的 调节 y& olm 两 部 分 。 在 可 解困 
的 调节 中 包含 调节 变量 Z 直接 的 调节 和 通过 中 介 的 调节 变量 M 间接 的 调节 。 因 此 ， 不 论 是 


误差 方差 满足 齐 性 的 条 件 还 是 误差 方差 依赖 于 X 的 情况 ， 我 们 都 给 出 了 测量 有 中 介 调 节 效 
应 的 效应 量 指标 。 公式 (13) - C15) 同时 适用 于 meMO 和 2meMO 模型 ,而 公式 (10) - (12) 


仅 适 用 于 2meMO 模型 。 


由 于 本 研究 拓展 的 2meMO 模 型 随机 部 分 参数 更 复杂 ， 极 大 似 然 估计 方法 会 遇 到 基本 假 
设 (误差 方差 齐 性 ) 被 违背 的 问题 。 考 虑 到 贝 叶 斯 估计 在 处 理 复杂 模型 ， 尤 其 是 含有 更 复杂 


的 随机 变异 模型 时 更 加 灵活 的 优势 (Wang & Preacher, 2015; Muthén & Asparouhov, 2012; 
Asparouhov & Muthén, 2020)， 我 们 使 用 贝 叶 斯 方法 估计 2meMO 中 的 参数 及 效应 量 。 为 了 便 


于 比较 ， 贝 叶 斯 估计 方法 也 用 于 meMO 模 型 的 参数 估计 。 在 指定 模型 和 先 验 分 布 后 ， 采 用 马 


尔 可 夫 链 蒙特 卡 洛 (Markov Chain Monte Carlo, MCMC ) Gibbs 采 样 获得 模型 参数 的 后 验 分 布 


(Gilks et al., 1996)。 关 于 Gibbs 采 样 收敛 性 的 检查 ， 使 用 Gelman-Rubin 势 能 缩减 (potential 


scale reduction, PSR) 统 计量 。 根 据 Brooks 和 Gelman(1998) 的 研究 , 如 果 PSR 统 计量 小 于 1 
则 认为 参数 估计 达到 收敛 。 对 于 收 和 敛 样本 , 我 们 将 参数 的 贝 叶 斯 估计 定义 为 其 后 验 分 布 


.05， 
的 均 


值 〈 称 为 后 验 均值 ， 实 际 中 也 可 用 中 位 数 ) ， 并 且 使 用 收敛 样本 的 标准 差 获得 后 验 标 准 


Bo 


(Song & Lee, 2012) 。 


4 模拟 研究 


基于 每 个 参数 后 验 分 布 2.5% 和 97.5% 的 分 位 数 ， 得 到 95% 可 信 区 间 (Credibility Interval, CI) 


本 研究 通过 模拟 研究 的 方法 , 比较 不 同 条 件 下 , meMO 和 2meMO 在 参数 估计 精度 和 效应 


四 | 


2meMO 模 型 的 必要 性 。 


量 估计 方面 的 表现 。 通 过 meMO 和 2meMO 的 比较 ， 验 证 方差 齐 性 条 件 不 满足 时 ， 拓 展 的 


4.1 模拟 设计 


对 于 方程 (2) - (5$)， 参 考 Wen 等 人 (2010) 以 及 Liu 等 人 (2019) 的 模拟 条 件 ， 参 数 真 值 的 设 


定 如 下 : C1) 调节 所 FRB ra 的 大 小 取 3 个 水 平 : 


0、0.2 和 0.4; (2) 随机 误差 ui 的 方差 o2。 


取 3 个 水 平 : 0、0.25 和 0.5; (3) 样本 量 取 4 个 水 平 : 100、200、500 和 1000。 注意 条 件 o2。 = 0 


时 ， 表 示 了 的 误差 方差 满足 齐 性 假设 ， 这 时 对 应 的 真 模型 为 neMO 模 型 。 其 


他 条 件 的 设 定 如 


TF: 将 截 距 yy。 和 ygo 设 置 为 0; HX >Y, Z>Y, M > YMZ > M 的 主 效应 均 设 置 为 0.4， 即 


Va1 = Yaz = Yeo = Ym1 = 0.4; 2 对 X 一 了 路 径 的 直接 调节 效应 设置 为 0.2， Bycs = 0.2。 自 变 


500 次 。 基 于 2meMO 模 型 ， 使 用 R 软 件 (R Core Team, 2016) 模拟 生成 数据 。 可 


C15) 可 以 得 到 对 应 于 yecz 和 ocxc 组 合 下 ， 


量 X， 调 节 变 量 Z2， 误 差 项 eyy; 和 ey; 服 从 标准 


效应 量 的 真 值 ， 结 果 见 表 1。 


表 1 模拟 设计 条 件 对 应 的 效应 量 真 值 


正 态 分 布 M0,1)。 共 有 36 个 条 件 ， 每 个 条 件 重复 


民 据 公式 (10) - 


Ycı 0 0.2 0.4 

o2, 0 0.25 0.5 0 0.25 05 0 0.25 05 

as 1 1 1 0.6622 0.6622 0.6622 0.4475 0.4475 0.4475 

Ve i 1 1 1 0.3378 0.3378 0.3378 0.1381 0.1381 0.1381 
PLP uo 0 0 0 0.0541 0.0541 0.0541 0.0884 0.0884 0.0884 
mo tot 1 0.1379 0.7410 0.6622 0.2128 0.1268 0.4475 0.2402 0.1641 
mo dir 1 0.1379 0.7410 0.3378 0.1086 0.0647 0.1381 0.0741 0.0507 
Pmemo 0 0 0 0.0541 0.0174 0.0103 0.0884 0.0474 0.0324 


4.2 数据 分 析 和 评价 指标 


对 于 2meMO 模 型 和 meMO 模 型 ， 使 用 JAGS(Plummer, 2015) 和 R 软 件 (R Core Team, 2016) 


中 的 程序 包 “R2jags” 进 行 贝 叶 斯 估计 ,共和 夫 代 15000 次 , 采用 最 后 5000 次 迭代 获得 参数 的 后 


验 分 布 , 即 内 燃 (burn in) 的 值 设 定 为 10000。 对 于 模型 中 参数 先 验 分 布 的 选择 , 基于 文献 建议 


(e.g., Song & Lee, 2012; Wang & Preacher, 


(Ymo Yui» Yao Var» Yaz Voor Ver Ye2) a 月 


H 4k 立正 态 分 布 的 先 验 ， 
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2015) ， 


x F 


H 系数 y= 
对 于 误差 方差 oz = 


(a2m,027,026) 采 用 独立 的 逆 伽 马 分 布 的 先 验 , 这 些 先 验 是 回归 分 析 和 结构 方程 模型 常用 的 半 
He 46 4e(Gelman et al.,2004)。 有 具体 而 言 ， 对 于 回归 系数 ， 设 定 先 验 分 布 为 : yx~N(0,109); 
对 于 误差 方差 , 设 定 先 验 分 布 为 : of~1Gamma(0.001,0.001)。 链 条 数 设 为 3, 为 了 降低 后 验 
分 布 间 的 自 相关 和 节约 存储 空间 ， 将 抽取 后 验 分 布 的 间隔 设 定 为 10〈 即 thin=10)。 感 兴趣 的 
参数 包括 中 介 调 节 效应 ymewo 和 中 介 调节 效应 的 差异 dijrewo， 取 2 = 1 来 评估 的 简单 meMO 
效应 (ymemo) 的 估计 ， 取 2 = 1 和 2 = 一 1 来 检验 调节 变量 Z 取 不 同 值 时 两 个 简单 meMO 效 应 
的 差异 (difnemo)。 

对 于 参数 yi,emo 和 difinemo， 使 用 参数 估计 的 偏差 (Bias)、 均 方 误差 (MSE)、95% 可 信 
区 间 的 覆盖 率 、 检 验 力 和 第 一 类 错误 5 个 指标 评估 参数 估计 的 精度 。 对 于 效应 量 , 采用 前 3 个 


指标 评估 其 估计 精度 。 同 时 采用 meMO 和 2meMO 模 型 进行 估计 ， 对 于 2meMO 模 型 ， 可 以 得 


BG? ror OD ar P2Quor Puos $dir 和 memo 的 估计 结果 ， 对 于 meMO 模 型 ， 可 
得 到 前 3 个 效应 量 的 估计 值 。 
偏差 (Bias) 和 均 方 误差 (MSE) 的 定义 如 下 : 
Bias(6) = 6 — 6, (16) 


MSE(6) = :Co (17) 


其 中 ，5 = YR Â RRI R=-500 次 重复 参数 估计 的 均值 ，B 为 第 "次 重复 的 参数 估计 值 。 

95% 可 信 区 间 的 覆盖 率 是 指 R=500 次 重复 中 贝 叶 斯 估计 得 到 的 后 验 分 布 中 第 2.5% 和 第 
97.5% 可 信 区 间 内 包含 参数 真 值 的 次 数 所 占 的 比例 。 检 验 力 是 在 参数 对 应 的 真 值 不 为 0 的 条 
件 下 ，R=500 次 重复 中 能 够 检验 出 参数 显著 不 为 零 的 次 数 所 占 的 比例 。 第 一 类 错误 是 在 参数 
对 应 的 真 值 为 0 的 条 件 下 ，R=500 次 重复 中 错误 的 检验 出 参数 显著 不 为 0 的 次 数 所 占 的 比例 。 


4. 3 结果 
对 于 每 一 组 模拟 条 件 , 模型 中 所 有 参数 对 应 的 PSR 统计 量 均 小 于 1.05, 说 明 对 于 meMO 
和 2meMO 模型 ， 贝 叶 斯 佑 计 均 不 存在 不 收敛 的 问题 。 


4. 3.1 ymemo 和 difmemo 的 估计 精度 


表 2 给 出 了 有 中 介 调 节 效 应 估计 值 ?emo 的 偏差 , 均 方 误差 , 95% 可 信 区 间 履 盖 率 和 拒绝 
率 。 对 于 拒绝 率 , Ymemo = 0 时 拒绝 率 的 值 代 表 第 一 类 错误 率 , Vmemo + 0 时 拒绝 率 的 值 代表 
检验 力 。 可 以 看 出 , 对 于 参数 估计 的 偏差 和 均 方 误差 两 种 模型 得 到 的 参数 估计 结果 几乎 没有 
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差异 〈 差 异 在 小 数 点 后 第 四 位 ) 。 在 满足 方差 齐 性 的 条 件 下 ，meMO 模 型 与 2meMO 模 型 得 到 
的 参数 估计 结果 也 基本 一 致 。 但 是 当 方差 齐 性 条 件 不 满足 时 ， 随 着 oz 方差 增加 ， 对 于 95% 可 
信 区 间 履 新 率 和 第 一 类 错误 的 结果 , 2meMO 模 型 的 结果 明显 更 有 优势 .对 于 检验 力 的 结果 ， 
meMO 在 小 样本 时 略 高 于 2meMO, 其 主要 原因 可 能 是 因为 这 些 条 件 下 , meMO 模 型 过 高 的 第 
一 类 错误 导致 检验 力 的 虚假 提高 。 随 着 样本 量 和 调节 回归 系数 y。 的 增加 ， 检 验 力 提高 。 
difineM0 与 ymeMo 所 得 结论 一 致 。 传 统 meMO 和 2meMO 模 型 的 偏差 都 很 小 。 在 方差 齐 性 
的 条 件 下 , 两 种 模型 中 Bias、MSE、 对 应 的 95% 可 信 区 间 的 覆盖 率 和 第 一 类 错误 的 差异 很 小 。 


然而 ， 在 误差 方差 非 齐 性 条 件 下 ，2meMO 对 difiemo 进 行 估计 更 加 准确 ， 同 


时 2meMO 对 应 


的 95% 可 信 区 间 的 履 盖 率 更 接近 95%。2meMO 方 法 在 控制 第 一 类 错误 率 方面 表现 更 好 ， 随 着 


03. 值 的 增加 ，meMO 的 第 一 类 错误 率 有 所 增 大 。 


4. 3. 2 效应 量 的 估计 结果 


对 于 总 体 中 存在 中 介 调 节 效 应 (ycs = 0.2 或 0.4) 时 ， 可 以 使 用 meMO 模 


型 和 2meMO 模 


型 ， 利 用 公式 (13) - (15) BBO or OO or 和 yo2wo 的 估计 值 ; 当 回归 系数 ci 对 


应 的 误差 方差 c2. + OM, 可 以 采用 2meMO 模 型 , 利用 公式 (10) - (12) 得 到 效应 量 @ wo tot 


pmo_air 和 gmemo 的 估计 值 。 效应 量 的 结果 表明 , 不 同 测量 指标 的 结果 表现 出 相同 的 模式 ， 所 


以 文中 仅 给 出 pWwyo 和 gpmemo 的 结果 。 


表 3 给 出 了 两 种 模型 得 到 的 p50 的 估计 值 。 样 本 量 越 大 ，o2, 的 值 越 小 ， 


Oana ttt 


HERE. Eog = 0 的 情况 下 ，meMO 模 型 和 2meMO 模 型 间 的 差异 很 小 。 但 是 ， 当 ax #0 


时 ，2meMO 的 结果 更 加 精确 ， 并 且 随 着 oxe 值 的 增加 ，2meMO 模 型 相对 于 meMO 模 型 的 优势 


变 得 越 来 越 明显 。 
表 4 给 出 了 2meMO 模型 得 到 的 效应 量 Bpmemo 的 估计 值 。 从 结果 可 以 看 


HH， 在 所 有 条 


IEF, MSE 值 都 很 小 ，pyemo 对 应 的 95% 可 信 区 间 的 履 盖 率 也 接近 95%。 随 着 oz2。 和 样本 


量 的 增加 ，gpmemo 估 计 的 精度 越 来 越 高 。 但 是 ， 值 得 注意 的 是 ， 当 样本 量 较 小 时 ，@,emo 
估计 的 偏差 较 大 ， 尤 其 是 对 于 较 小 的 有 中 介 的 调节 效应 而 言 。 因 此 ， 在 实践 中 ， 为 了 获得 


对 pmemo 更 准确 的 估计 ， 可 能 需要 更 大 的 样本 量 ， 或 者 类 似 于 传统 回归 分 析 中 对 R Bf 


计 ， 需 要 对 其 偏 度 进行 校正 。 
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表 2 两 种 模型 得 到 的 ymewo 的 偏差 ， 均 方 误差 ，95% 可 信 区 间 履 盖 率 和 拒绝 率 


meMO 2meMO 


2 
o, 真 值 “一 一 
ma 人 100 200 500 1000 100 200 500 1000 


Bias 0 0 0 0.002 0.001 0.000 -0.001 0.002 0.001 0.000 -0.001 


0.25 0 -0.002 -0.001 0.000 0.001 -0.001 -0.001 0.000 0.001 
0.5 0 0.003 0.000 0.001 -0.001 0.003 0.000 0.001 -0.001 
0.2 0 0.08 -0.002 -0.002 0.000 0.001 -0.002 -0.002 0.000 0.001 


0.25 0.08 -0.006 0.000 0.000 0.000 -0.006 0.000 0.000 0.000 
0.5 0.08 -0.002 0.002 -0.002 0.001 -0.001 0.000 -0.002 0.002 
0.4 0 0.10 0.000 0.001 0.001 0.000 0.000 0.001 0.001 0.000 
0.25 0.16 -0.004 -0.001 0.002 0.000 -0.004 -0.001 0.002 0.000 
0.5 0.16 -0.002 -0.004 0.000 0.002 -0.002 -0.003 0.001 0.002 


MSE 0 0 0 0.002 0.001 0.000 0.000 0.002 0.001 0.000 0.000 
0.25 0 0.003 0.001 0.001 0.000 0.003 0.001 0.001 0.000 

0.5 0 0.004 0.002 0.001 0.000 0.004 0.002 0.001 0.000 

0.2 0 0.08 0.002 0.001 0.000 0.000 0.002 0.001 0.000 0.000 


0.25 0.08 0.004 0.002 0.001 0.000 0.003 0.002 0.001 0.000 
0.5 0.08 0.005 0.002 0.001 0.001 0.005 0.002 0.001 0.000 
0.4 0 0.16 0.004 0.002 0.001 0.000 0.004 0.002 0.001 0.000 
0.25 0.16 0.005 0.002 0.001 0.000 0.005 0.002 0.001 0.000 
0.5 0.16 0.006 0.003 0.001 0.001 0.005 0.003 0.001 0.001 


ime 0 0 0 0.952 0.942 0.942 0.942 0.952 0.956 0.944 0.940 
0.25 0 0.908 0.920 0.918 0.872 0.940 0.956 0.942 0.932 

0.5 0 0.912 0.874 0.852 0.870 0.948 0.952 0.936 0.948 

02 0 0.08 0.962 0.950 0.966 0.924 0.968 0.954 0.972 0.926 


0.25 0.08 0.920 0.918 0.922 0.908 0.930 0.944 0.956 0.946 
0.5 0.08 0.894 0.882 0.884 0.882 0.928 0.946 0.936 0.944 
0.4 0 0.16 0.960 0.948 0.958 0.950 0.962 0.950 0.958 0.946 
0.25 0.16 0.935 0.937 0.930 0.933 0.960 0.948 0.942 0.962 
0.5 0.16 0.926 0.928 0.896 0.882 0.946 0.954 0.940 0.940 


拒绝 率 


© 
© 
© 


0.048 0.058 0.058 0.058 0.048 0.044 0.056 0.060 

025 0 0.092 0.080 0.082 0.129 0.060 0.044 0.058 0.068 

0.5 0 0.088 0.126 0.148 0.130 0.052 0.048 0.064 0.052 

02 0 0.08 0.438 0.762 0.998 1.000 0.399 0.760 0.996 1.000 

0.25 0.08 0.390 0.650 0.970 1.000 0.314 0.606 0.950 1.000 

0.5 0.08 0.362 0.600 0.882 0.994 0.269 0.488 0.838 0.990 

04 0 0.16 0.938 1.000 1.000 1.000 0.926 1.000 1.000 1.000 

0.25 0.16 0.834 0.996 1.000 1.000 0.808 0.952 1.000 1.000 

0.5 0.16 0.794 0.970 1.000 1.000 0.730 0.964 1.000 1.000 

TE: 覆盖 率 低 于 0.9， 第 一 类 错误 率 高 于 8% 的 值 用 粗 体 字 表示 。ymemo = 0 时 拒绝 率 的 值 代 表 第 一 类 错误 
率 ，ymemo + 0 时 拒绝 率 的 值 代表 检验 力 。 
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表 3 两 种 模型 得 到 的 效应 量 pYjyo 估 计 的 偏差 ， 均 方 误差 ， 95% 可 信 区 间 覆 盖 率 


> moME 2moME 
Neil Cuc 真 值 
100 200 500 1000 100 200 500 1000 
Bias 02 0 0.054 0.014 0.006 0.002 0.002 0.015 0.006 0.003 0.002 


0.25 0.054 0.015 0.009 0.004 0.001 0.016 0.010 0.004 0.001 
0.5 0.054 0.023 0.013 0.003 0.003 0.022 0.012 0.003 0.003 


04 0 0.088 0.003 0.001 0.001 0.001 0.004 0.001 0.001 0.001 
0.25 0.088 0.000 0.001 0.002 -0.001 0.000 0.001 0.002 -0.001 
0.5 0.088 0.001 -0.002 0.002 0.001 0.002 -0.002 0.001 0.001 


MSE 02 0 0.054 0.003 0.001 0.001 0.000 0.003 0.001 0.001 0.000 
0.25 0.054 0.003 0.002 0.001 0.000 0.003 0.002 0.001 0.000 
0.5 0.054 0.004 0.002 0.001 0.001 0.003 0.002 0.001 0.000 


04 0 0.088 0.002 0.001 0.000 0.000 0.002 0.001 0.000 0.000 
0.25 0.088 0.003 0.001 0.001 0.000 0.003 0.001 0.001 0.000 
0.5 0.088 0.003 0.002 0.001 0.000 0.003 0.001 0.001 0.000 


Pez 


aK 02 0 0.054 0.962 0.940 0.956 0.926 0.958 0.944 0.954 0.922 
0.25 0.054 0.958 0.926 0.922 0.914 0.970 0.950 0.962 0.948 
0.5 0.054 0.958 0.902 0.882 0.874 0.980 0.954 0.952 0.934 


04 0 0.088 0.968 0.952 0.948 0.954 0.966 0.952 0.948 0.950 
0.25 0.088 0.964 0.944 0.948 0.952 0.964 0.944 0.948 0.952 
0.5 0.088 0.932 0.908 0.914 0.874 0.950 0.940 0.956 0.932 


Es 覆盖 率 低 于 0.9 的 值 用 粗 体 字 表 示 。 
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表 4 2meMO 模型 得 到 的 效应 量 由 exwo 估 计 的 偏差 ， 均 方 误差 ，95% 可 信 区 间 履 盖 率 


Yc1 aže HE 100 200 500 1000 

Bias 0 0 0 0.040 0.024 0.013 0.007 
0.25 0 0.030 0.016 0.005 0.002 

0.5 0 0.023 0.009 0.003 0.001 

0.2 0 0.0541 -0.002 -0.005 -0.006 -0.004 

0.25 0.0174 0.021 0.016 0.006 0.002 

0.5 0.0104 0.020 0.010 0.002 0.002 

0.4 0 0.0884 -0.009 -0.008 -0.005 -0.004 

0.25 0.0474 0.013 0.011 0.006 0.002 

0.5 0.0324 0.017 0.007 0.003 0.002 

MSE 0 0 0 0.003 0.001 0.000 0.000 
0.25 0 0.002 0.001 0.000 0.000 

0.5 0 0.001 0.000 0.000 0.000 

0.2 0 0.0541 0.002 0.001 0.000 0.000 

0.25 0.0174 0.001 0.001 0.000 0.000 

0.5 0.0104 0.001 0.001 0.000 0.000 

0.4 0 0.0884 0.002 0.001 0.000 0.000 

0.25 0.0474 0.002 0.001 0.000 0.000 

0.5 0.0324 0.002 0.001 0.000 0.000 

imi 0.2 0 0.0541 0.932 0.936 0.948 0.922 
0.25 0.0174 0.978 0.954 0.954 0.936 

0.5 0.0104 0.976 0.964 0.940 0.948 

0.4 0 0.0884 0.946 0.946 0.950 0.932 

0.25 0.0474 0.962 0.946 0.936 0.936 

0.5 0.0324 0.960 0.962 0.932 0.948 


VE: 


对 于 2meMO 模型 ， 


于 方差 估计 大 于 0， 置信 区 间 的 覆盖 率 均 不 包含 0。 


5 应 用 案例 


5.1 数据 和 模型 


所 


为 了 便于 2meMO 模 型 与 meMO 模 型 的 比较 ， 这 里 我 们 采用 Kwan 和 Chan (2018) & 


究 中 


的 数据 (下载 地 址 :http://dx.doi.org/10.1037/met0000160.supp) 。 其 数据 来 源 于 2012 年 


OECD (Organization for Economic Cooperation and Development) 组 织 的 PISA 测 试 (the Program 


for International Student Assessment, PISA), 


该 数据 


能 (MEFF)， 调 节 变 量 社 会 经 济 地 


本 量 为 3047， 分 析 的 变量 包括 : Wl 


YESCS) 以 及 中 介 的 调节 变量 家 庭 教育 资源 (HRES)。 模 型 


具体 介绍 见 Kwan 和 Chan (2018) 研 究 中 的 介绍 


o 


变量 数学 学 业 成 绩 (MATH)， 自 变量 数学 自我 效 


假设 : ESCS 可 以 调节 MEFF 与 MATH 之 间 的 关系 ，ESCS 较 低 的 学 生 ，MEFF 对 MATH 的 影响 


较 强 ， 而 ESCS 较 高 的 学 生 ，MEFF 对 MATH 的 
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i= 
We 


响 较 弱 ; 进一步 ，ESCS 对 MEFF 与 MATH 之 


间 关 系 的 调节 作用 被 HRES 所 中 介 ，ESCS 高 的 学 生 拥 有 较 多 的 HRES,， 所 以 MEFF 对 MATH 的 
影响 较 弱 。 对 应 的 概念 模型 见 图 4， 分 析 的 模型 包括 传统 meMO 模 型 和 2meMO 模 型 ， 分 析 过 
程 中 对 MEFF、ESCS 和 HRES 做 均值 中 心 化 处 理 。 为 了 帮助 应 用 者 理解 和 便于 应 用 者 使 用 分 
析 程 序 ， 我 们 在 附录 中 给 出 了 默认 先 验 设置 时 Mplus8.3 的 语句 (附录 1) 。 需 要 再 次 强调 的 是 
这 里 的 数据 是 单 层次 结构 , 在 Mplus 中 我 们 借助 于 多 层 模型 采用 贝 叶 斯 估计 方法 实现 2meMO 


的 参数 估计 ， 一 个 个 体 就 是 一 个 cluster。 


MEFF 


4 家 庭 教育 资源 对 社会 经 济 地 位 与 数学 自我 效能 与 数学 成 绩 之 间 中 介 调 节 效 应 模型 


5.2 分 析 结 果 


对 于 meMO 模 型 和 2meMO 模 型 ， 在 默认 先 验 设置 条 件 下 模型 中 所 有 参数 对 应 的 PSR 统 计 
量 均 介 于 1.00 到 1.02 之 间 , 进一步 检查 3 条 不 同 链 的 轨迹 图 发 现 不 同 链条 收敛 到 同样 的 后 验 分 
布 。 表 6 给 出 了 两 个 模型 的 参数 估计 结果 、 中 介 的 调节 效应 估计 值 和 本 研究 提出 的 效应 量 指 


标 。2meMO 模 型 对 应 的 DIC (Deviance Information Criteria) 为 78000.4，meMO 模 型 对 应 的 


DIC 为 78069.2, 两 个 模型 对 应 的 DIC 的 差异 和 DIC=68.8, 2meMO 模 型 的 DIC 较 小 , 根据 Cain 和 


Zhang(2019) 所 建议 的 ADIC>7 的 模型 比较 标准 ，2meMO 模 型 的 拟 合 优 于 meMO 模 型 。 
从 表 5 的 结果 可 以 看 出 , 两 个 模型 得 到 的 参数 估计 结果 在 大 多 数 参数 上 相差 不 大 .meMO 


和 2meMO 估 计 得 到 的 有 中 介 调 节 效 应 的 后 验 均 值 分 别 为 -2.019 和 -2.165， 对 应 的 95%CI 分 别 


为 (-3.439, -0.614) 和 (-3.757, -0.673)，CI 不 包含 0， 说 明 有 中 介 的 调节 效应 显著 。2meMO 模 型 
结果 表明 , 第 一 层 数学 自我 效能 对 数学 成 绩 影响 的 总 变异 的 后 验 均值 var(ci) = 484.239, 95% 
的 CI 为 (229.158, 737.546); 数学 自我 效能 对 数学 成 绩 影响 的 变异 中 可 以 被 调节 变量 社会 经 济 
地 位 和 中 介 的 调节 变量 家 庭 教育 资源 联合 解释 部 分 (ycz 十 Yerym1)>Var(Zi) + yê oem HY a om 


均值 为 46.643，95% 的 CI 为 (13.412, 83.211)。 由 贝 叶 斯 估计 得 到 的 效应 量 @Bjyo tot Pmo_air Fl 


mnemo 的 后 验 均 值 分 别 为 0.059,0.025 和 0.010，95% CI 分 别 为 (0.007, 0.122), (0.000, 0.068) 和 
(0.000, 0.023)， 说 明 在 数学 成 绩 对 数学 自我 效能 影响 的 总 变异 中 ， 社 会 经 济 地 位 调节 效应 的 
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总 解释 率 的 后 验 均 值 为 5.9%, 社会 经 济 地 位 直接 调节 效应 的 解释 率 的 后 验 均值 为 2.5%, 间接 
调节 效应 〈 即 有 中 介 的 调节 效应 ) 的 解释 率 的 后 验 均值 为 1.0%。 分 母 上 去 掉 误差 uci 的 方差 


后 ， 得 到 效应 量 6Ve 


MO_tot’ 


PO) 


R jw 和 BY2io 的 后 验 均值 的 估计 值 分 别 为 0.637，0.275 和 0.106， 


95% CI 分 别 为 (0.296,1.000)，(0.000,0.667) 和 (0.000,0.206)， 说 明 在 不 考虑 数学 自我 效能 一 数 


学 学 业 成 绩 影响 系数 的 随机 误差 (wci) IY, FEA 


会 经 济 地 位 和 家 庭 教 


育 资源 联合 解释 的 数学 


自我 效能 对 数学 成 绩 影响 的 变异 中 ， 社 会 经 济 地 位 总 解释 率 的 后 验 均值 为 63.7%， 社 会 经 济 
地 位 直接 调节 效应 的 解释 率 的 后 验 均值 为 27.5%， 间 接 调节 效应 的 解释 率 的 后 验 均值 为 


10.6%. ASF! 


也 从 另 一 角度 说 明了 2meMO 模 型 ， 


, meno SO? ot SUA LE RRA, HERR 


考虑 X 一 了 随机 误差 的 意义 。 


表 5 meMO 和 2meMO 模 型 参数 估计 结果 


因为 所 用 分 母 的 差异 ， 


统计 量 meMO 2meMO 
SUV Ee ea au 
后 验 均值 ” 后 验 标准 差 MJE ” 后 验 标准 差 
路 径 系数 
ESCS—HRES (fm1) 0.484*** 0.011 0.484*** 0.011 
MEFF—>MATH (7.0) 46.438*** 1,385 48.436*** 1,571 
ESCS—MATH (fa2) 14.621*** 1.687 14.368*** 1.676 
HRES—MATH (741) 3.492* 1.672 3.418* 1.644 
ESCSxMEFF—>MATH (7.2) -2.911* 1.501 -3.284* 1.640 
HRESxMEFF—>MATH (7.1) -4.171** 1.478 -4.471** 1.631 
ESCSxMEFF—>HRESxMEFF (11) 0.484*** 0.011 0.484*** 0.011 
LET 6197.923 160.921 5668.715 189.195 
a 0.749 0.019 0.749 0.019 
Te - - 484.239 130.388 
有 中 介 调 节 meMO 
Pume -2.019 0.718 -2.165 0.791 
效应 量 
Po tot = 0.059 0.032 
bmo_air - - 0.025 0.022 
memo - - 0.010 0.007 
BO sot 0.623 0.193 0.637 0.192 


z0) 
PR air 0.261 0 


pi) 
EL mo 0.110 0 


.203 0.275 


057 0.106 


本 例 中 ， 通 过 模型 拟 合 的 比较 我 们 倾向 于 采用 2meMO， 研 究 者 可 以 基于 感 兴趣 的 研究 


问题 ， 在 公式 〈10) - (15) 所 定义 的 6 种 不 同 的 效应 量 中 选择 合适 的 效应 量 ， 解 释 调节 变量 
对 X 一 Y 之 间 关 系 变异 的 解释 。 但 如 果实 际 中 ， 我 们 倾向 于 选择 meMO 模 型 ， 则 根据 感 兴趣 


的 研究 问题 ， 在 公式 13) - C15) 中 选择 合适 的 效应 量 


型 的 转换 过 程 中 ， 我 们 强调 模型 中 隐 含 了 的 对 X > XM 


。 另 外 ， 虽 然 从 概念 模型 到 统计 模 
和 XZ 一 XM 两 条 路 径 系 数 的 限定 ， 但 


在 实际 中 ， 这 两 个 限定 条 件 对 于 参数 估计 结果 的 影响 很 小 ， 甚 至 可 以 忽略 。 因 此 ,为 了 简化 


操作 ， 正 如 Kwan 和 Chan(2018) 指出 的 实际 中 也 可 以 久 


6 讨论 与 建议 


本 文 基于 单 层 的 数据 拓展 了 meMO 模型 ， 采 用 两 层 模 型 的 


略 这 两 个 限定 。 


EE 架 定 义 了 可 处 理 误差 方差 


非 齐 性 的 2meMO 模型 , 基于 此 模型 将 对 了 影响 系数 变异 的 来 源 进行 分 解 , 并 由 此 定义 了 


可 以 测量 总 调节 效应 、 直 接 调 节 效 应 和 有 中 介 的 调节 效应 大 小 的 效应 量 指标 ,。 进一步 将 贝 叶 
斯 方法 用 于 2meMO 模型 的 参数 估计 ,并 通过 Monte Carlo 模拟 研究 考查 了 参数 估计 的 精度 ， 
并 与 传统 的 meMO 模型 进行 了 比较 ,验证 了 2meMO 模型 在 处 理 方差 非 齐 性 条 件 下 的 优势 。 

根据 模拟 研究 的 结果 , 在 方差 满足 齐 性 的 假设 时 , meMO 和 2meMO 得 到 的 结果 几乎 没 


有 差异 ， 也 就 是 说 即便 在 传统 meMO 的 假设 下 , 2meMO 也 能 得 到 无 偏 的 精度 较 高 的 估计 结 


果 。 但 是 在 方差 齐 性 的 条 件 不 满足 时 ，2meMO 得 到 的 参数 估计 结果 在 95% 可 信 区 间 覆 盖 率 
和 控制 第 一 类 错误 上 的 表现 优 于 meMO 模型 ， 这 一 结果 与 方差 非 齐 性 的 潜 变 量 交互 作用 模 
型 (Liu et al., 2019) 和 观测 变量 的 调节 效应 回归 分 析 模 型 得 到 的 结果 一 致 (Yang & Yuan, 2016). 
2meMO 模型 的 拓展 为 调节 效应 模型 中 常见 的 方差 齐 性 假设 条 件 的 违背 (Aguinis & Pierce, 


1998; Alexander & DeShon, 1994) 提 供 了 解决 思路 。 基 于 贝 叶 斯 估计 方法 的 可 实现 性 以 及 模 


拟 研究 的 结果 ， 我 们 推荐 研究 者 使 用 2meMO 模型 处 理 方差 非 齐 性 的 情况 ， 


感 兴趣 的 研究 问题 相对 应 的 效应 量 指标 。 


在 研究 中 报告 与 


在 有 中 介 的 调节 效应 的 分 析 中 , 本 文 提出 的 2meMO 模型 和 对 应 的 效应 量 指标 有 以 下 几 
方面 的 优点 : (1) 采用 两 层 模型 定义 meMO 统计 模型 的 方式 ， 在 概念 模型 的 理解 和 效应 量 


的 定义 上 为 如 何 区 分 变量 的 角色 提供 了 更 为 一 般 化 的 术 


EZR. FY DUKE RAY A 
量 关 系 的 变化 ， 定 义 为 层 1 的 模型 ,然后 根据 调节 变量 的 定义 ,在 层 2 模型 中 解释 调节 变量 
18 


的 自 变量 和 因 变 


对 因 变 量 和 自 变 量 之 间 关 系 的 影响 ， 以 及 调节 变量 和 中 介 的 调节 变量 之 间 的 关系 。 这 样 分 层 
次 的 模型 表达 , 合理 且 有 效 地 区 分 了 自 变量 、 因 变量 、 调 节 变 量 和 中 介 的 调节 变量 的 不 同 角 
色 。(2) 文中 提出 的 效应 量 指标 是 基于 自 变量 和 因 变 量 之 间 总 变异 〈 即 var(c)) 的 分 解 所 得 
到 的 ， 这 就 为 紧密 围绕 “与 了 之 间 关 系 的 调节 ”这 一 关键 问题 ， 定 义 合适 的 效应 量 提供 了 
基础 ， 基 于 此 所 定义 的 效应 量 是 围绕 meMO 分 析 所 感 兴趣 的 问题 提出 的 ， 更 具 可 解释 性 ， 
同时 也 更 符合 效应 量 概念 的 核心 内 涵 (e.g., Kelley & Preacher, 2012). (3) 基于 2meMO 提出 
的 效应 量 指标 ， 不 仅 适 用 于 方差 齐 性 而 且 适 用 于 方差 非 齐 性 的 条 件 ， 有 更 广 的 适用 范围 。 研 
究 者 在 使 用 过 程 中 不 用 担心 由 于 实际 数据 不 满足 方差 齐 性 而 带 来 的 估计 偏差 , 检验 力 降 低 等 
问题 。(4) 2meMO 模型 的 参数 估计 很 容易 在 Gibbs 取样 的 框架 下 ， 借 用 已 有 的 免费 贝 叶 斯 
估计 程序 实现 。 如 2meMO 的 参数 估计 很 容易 通过 常用 的 JAGS (Plummer, 2015)、WinBugs 
(Lunn et al., 2012) 和 Mplus 等 软件 实现 。(S$) 更 重要 的 是 ， 本 文 所 提出 的 有 中 介 的 调节 模型 
的 效应 量 指标 可 满足 应 用 研究 者 报告 效应 量 的 需求 ， 弥 补 了 目前 尚 没有 合适 的 测量 meMO 


lAr 
14v1 


| 
y 


S 影响 的 效应 量 指标 的 局 限 。 

© 另外 ， 本 研究 定义 模型 的 方法 和 效应 量 的 思路 具有 拓展 性 。 本 文 提出 的 定义 2meMO 模 
= 型 的 思路 ， 可 以 在 概念 上 为 如 何 有 效 区 分 meMO 和 moME 两 类 模型 提供 一 些 借鉴 。 从 概念 模 
> 型 来 看 ，2meMO 模 型 和 moME-I 强 调 的 重点 不 同 ， 两 者 并 不 等 价 ，2meMO 模 型 是 为 了 解释 


| 


a!) 


《 
` 


调节 效应 是 如 何 通过 中 介 的 调节 变量 发 挥 作用 的 ; 而 moME-I 是 为 了 解释 中 介 效 应 是 如 何 依 
赖 于 调节 变量 取 值 的 。 然 而 ,传统 的 方差 齐 性 假设 下 由 于 两 者 对 应 的 统计 模型 相同 , 在 结 
Z 解释 上 很 难 有 效 帮 助 应 用 研究 者 区 分 两 类 模型 的 本 质 区 别 。 在 两 层 建 模 的 框架 下 ， 对 于 
moME 模 型 可 以 借鉴 本 文 的 建 模 思路 ， 在 层 1 定义 与 中 介 效应 有 关 的 路 径 Ca 和 2 ) ， 而 在 层 2 
解释 中 介 路 径 中 某 个 阶段 的 效应 被 调节 。 例 如 , 含有 第 二 阶段 被 调节 的 中 介 模 型 可 以 定义 为 : 
层 1: 


Mi = QiXi + emi 
Y; = di + b;M; + ¢;X; + & 
层 2: 
ai = Yao 
di = Yao + Var Zi + Uai 
bi = Yoo + Yor Zi + Upi 


Ci = Yeo 
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上 述 模 型 合并 后 的 表达 式 中 包含 Miupi; 这 一 项 ,而 2meMO 模 型 中 包含 的 是 Xiuc:， 这 样 统计 模 
型 的 定义 上 可 将 两 者 加 以 区 分 。 即 便 上 述 模型 中 不 包含 随机 误差 项 wy:， 这 一 两 层 模 型 建构 
的 思路 也 有 利于 研究 者 区 分 两 类 模型 解释 上 的 区 别 。 未 来 研究 中 可 以 借鉴 本 文 定 义 效应 量 的 
思路 , 基于 被 调节 的 中 介 所 关注 的 问题 发 展 合适 的 效应 量 指标 ,以 解决 有 中 介 的 调节 模型 与 
第 二 阶段 被 调节 的 中 介 模 型 经 常 出 现 的 结果 解释 上 被 混淆 的 问题 。 


Lachowicz, Preacher 和 Kelley (2018) 讨论 了 合理 的 效应 量 指标 应 该 满足 的 特性 (Kelley & 


Preacher, 2012; Preacher & Kelley, 2011; Wen & Fan, 2015)。 本 研究 所 定义 的 效应 量 指标 具有 
以 下 特性 。 (1) 效应 量 的 定义 与 neMO 的 概念 模型 匹配 ， 在 变异 分 解 的 框架 下 ， 效 应 量具 有 
可 解释 的 尺度 ， 可 以 回答 有 中 介 调 节 对 总 调节 变异 方差 的 贡献 率 。 (2) 效应 量 指标 不 依赖 


于 样本 量 ， 且 基于 贝 叶 斯 估计 的 后 验 分 布 可 得 到 效应 量 的 可 信 区 间 (CD 。 (3) 6? soe Fl 


Puo tot 的 取 值 介 于 0 和 1 之 间 。 大 多 数 情况 下 ，%go aire Puo Pmo.airMOmemo HHUA tH 


于 0 和 1 之 间 ; 但 在 实际 应 用 中 , 可 能 会 出 现 这 四 个 效应 量 的 数值 大 于 1 (Lutz, 1983; Maassen 
& Bakker, 2001; Tzelgov & Henik, 1991)， 此 时 说 明 Z 对 X 一 7 的 直接 调节 Ve) 和 通过 M 对 X 

一 7 的 间接 调节 (yerymi) 方向 相反 ， 出 现 了 变量 之 间 的 抑制 效应 Csuppression effect). 
虽然 所 提出 来 的 效应 量具 有 很 多 优点 ， 但 在 应 用 时 也 应 关注 其 局 限 性 。 首 先 ， 类 似 于 方 
差分 析 中 的 偏 *，var(cj) 是 定义 效应 量 的 基础 (分母)， 即 效应 量 的 解释 都 是 以 var(ci;) 作 为 比 
较 基 准 的 ， 因 此 这 些 效 应 量 本 质 上 是 一 个 相对 的 效应 量 指标 。 因此， 在 实际 应 用 中 ， 在 报告 
效应 量 的 时 候 应 同时 报告 分 母 的 大 小 ， 以 避免 过 度 硅 大 meMO 效 应 ， 同 时 也 有 助 于 了 解 在 X 
Fæ 对 7 的 影响 中 , var(ci) 的 大 小 是 否 有 实际 的 意义 。 研究 者 应 该 在 关注 var(ci) 本 身 大 小 的 前 提 下 ， 
进一步 解释 meMO 效 应 的 大 小 。 另 外 ， 我 们 建议 在 报告 效应 量 的 时 候 ， 同 时 报告 本 研究 提出 
的 多 个 效应 量 指标 ,为 全 面 了 解 变 量 的 影响 关系 和 强度 提供 更 多 的 补充 信息 。 同时， 应 该 注 

意 基于 变异 分 解 所 定义 的 直接 调节 效应 和 间接 调解 效应 的 效应 量 指标 不 具有 可 加 性 。 

大 量 的 研究 证 明 不 同 先 验 分 布 会 影响 参数 估计 的 精度 (Yuan & MacKinnon, 2009) , 因此 ， 
实际 中 如 何 选择 合适 的 先 验 分 布 是 研究 者 应 该 注意 的 问题 .本 研究 选取 了 最 常用 也 是 研究 者 
广泛 推荐 的 无 信息 先 验 (Browne & Draper, 2006) ， 结 果 证 明 即 使 是 无 信息 的 先 验 也 能 够 得 到 
准确 的 参数 估计 结果 , 但 是 这 一 结论 在 实际 应 用 中 并 不 一 定 具有 普 适 性 。 已 有 研究 表明 合适 
的 有 信息 先 验 可 以 使 得 参数 估计 精度 更 高 , Zondervan-Zwijnenburg 等 人 (2017) 提供 了 如 何 收 
集 先 验 信息 的 指导 。 如 果 研 究 者 可 以 得 到 合适 的 有 信息 的 先 验 , 我 们 也 鼓励 采用 有 信息 的 先 


验 以 提高 参数 估计 的 精度 。 然 而 ， 为 了 避免 由 于 特定 先 验 设置 带 来 的 主观 性 ,应 当 考虑 采用 
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介 


敏感 性 分 析 进 一 步 确认 贝 叶 斯 估计 结果 并 不 严重 依赖 于 某 个 特定 的 先 验 分 布 (van de Schoot 
et al., 2017; Zondervan-Zwijnenburg et al., 2017). 

本 研究 只 考虑 了 含有 一 个 调节 变量 Z 和 一 个 中 介 的 调节 变量 M 的 情况 ， 对 于 方差 非 齐 
性 也 只 考虑 了 结果 变量 了 的 方差 非 齐 性 ， 未 来 研究 中 可 以 考虑 更 为 复杂 的 meMO 模型 和 更 
复杂 的 方差 非 齐 性 的 情况 。 另 外 ， 对 于 本 研究 提出 的 效应 量 指标 ， 如 何 结合 效应 量 的 特点 ， 
给 出 效应 量 “ 小 入 “中兴 “ 大 ”的 判断 标准 ， 也 是 未 来 值得 进一步 研究 的 问题 。 最 后 ， 模 型 
中 被 忽略 的 变量 、 变 量 的 测量 误差 和 其 他 影响 回归 系数 估计 精度 的 因素 也 会 对 效应 量 的 估计 


精度 产生 影响 ， 未 来 研究 需要 对 这 些 影响 因素 进行 更 为 系统 和 深入 的 探讨 。 
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附录 1: meMO 和 2meMO 模型 贝 叶 斯 估计 的 Mplus 语句 
1. meMO 模型 
TITLE: mediated moderation model 
DATA: FILE = memo.csv; 
VARIABLE: NAMES ARE ID MATH ESCS HRES MEFF CESCS CHRES CMEFF XCESCS XCHRES; 
! CESCS CHRES CMEFF 为 均值 中 心 化 的 变量 ; 
! XCESCS= CMEFF*CESCS, XCHRES= CMEFF* CHRES; 
USEVARIABLES = MATH CESCS CHRES CMEFF XCESCS XCHRES; 
ANALYSIS: ESTIMATOR IS BAYES; POINT =MEAN; CHAINS =3; PROCESSORS =3; 
THIN = 10; BITERATIONS=(10000); BCONVERGENCE = 0.025; 
MODEL: 
MATH ON CMEFF (rc0); 
MATH ON CHRES (rdy1); 
MATH ON CESCS (rdy2); 
CHRES ON CESCS (rm1); 
MATH ON XCHRES (rc1); 
MATH ON XCESCS (rc2); 
XCHRES ON CMEFF (rm0); 
XCHRES ON XCESCS (rml); 
[CHRES MATH] (rm0 ry0); 
CHRES MATH (sig_em, sig_ey); 
CMEFF WITH CESCS; 
CMEFF WITH XCESCS; 
CESCS WITH XCESCS; 


MODEL CONSTRAINT: 
NEW(var_z vl v2 v3 v5 phi fl phi f2 phi _f3 ); 
var z=0.935; !z 的 方差 ; 
vl = (re2+re1*rm1)*(re2+1e1*rm1)*var_z; 
v2 =re2*re2*var_z; 
v3 =rel*rm1*rcl*rm1*var_z; 
v5 =vl4rel*rel*sig_em; 
phi fl =v1/5; 
phi f2 = v2/v5; 
phi f3 = v3/VS; 
MODEL CONSTRAINT: 
NEW(memo memol memo2 dif); 
memo=rcl*rml; 
memol =rcl*rml*0.967;  !z 的 标准 差 为 0.967; 


memo2 = rcl*rm1*(-0.967); 


dif = memol-memo2; 
PLOT: TYPE = PLOT2; 
OUTPUT: CINTERVAL(hpd) TECH8; 


25 


2. 2meMO 模型 的 Mplus 语句 
TITLE: two-level mediated moderation model with single-level data (2moME) 
DATA: FILE = memo.csv; 
VARIABLE: NAMES ARE ID MATH ESCS HRES MEFF CESCS CHRES CMEFF XCESCS XCHRES; 
USEVARIABLES = MATH CESCS CHRES CMEFF XCESCS XCHRES; 
CLUSTER = ID;。”! 对 于 单 层 数据 每 一 个 组 只 有 一 名 被 试 
WITHIN = MATH CESCS CHRES CMEFF XCESCS XCHRES; ! 单 层 数据 中 指定 所 有 变量 为 组 内 变量 
ANALYSIS: TYPE IS TWOLEVEL RANDOM; 
ESTIMATOR IS BAYES; 
POINT =MEAN; 
CHAINS =3; 
PROCESSORS =3; 
THIN = 10; 
BITERATIONS=(10000); 
BCONVERGENCE = 0.025; 


MODEL: 
%WITHIN% 

c | MATH ON CMEFF; 
MATH ON CHRES (rdy1); 
MATH ON CESCS (rdy2); 
CHRES ON CESCS (rm1); 
MATH ON XCHRES (rc1); 
MATH ON XCESCS (rc2); 
XCHRES ON CMEFF (rm0); 
XCHRES ON XCESCS (rml); 
[CHRES] (rm0); 

[MATH] (ry0); 

CHRES (sig_em); 
MATH; 

CMEFF WITH CESCS; 
CMEFF WITH XCESCS; 
CESCS WITH XCESCS; 


%BETWEEN% 
[c] (rc0); 
c (sig2_uc); 


MODEL CONSTRAINT: 
NEW(var_z vl v2 v3 v4 v5 phi_l phi_2 phi_3 phi_fl phi_f2 phi_f3 ); 


var_z =0.935; ! the variance of z; 

vl = (rc2+rc1*rm1)*(rc2+rc1*rm1)*var_z; 
v2 =rce2*re2*var_z; 

v3 =rel*rm1*rcel*rm1*var_z; 


v4 = vl4rel*rcl*sig_em+sig2_uc; 
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v5 = vl4rel*re1*sig_em; 
phi_1 = v1/v4; 
phi_2 = v2/v4; 
phi_3 = v3/v4; 
phi_fl = v1/v5; 
phi_f2 = v2/v5; 
phi_f3 = v3/v5; 


MODEL CONSTRAINT: 
NEW(memo memo! memo? dif); 
memo = rcl*rm1; 
memol = rcl*rm1*0.967; 
memo2 = rel *rm1*(-0.967); 
dif = memol-memo?; 

PLOT: TYPE = PLOT2; 

OUTPUT: 

CINTERVAL(hpd); 
TECH8; 
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Abstract 


Mediation and moderation analyses are commonly used methods for studying the relationship 
between an independent variable (X) and a dependent variable (Y) in conducting empirical research. 
To better understand the relationships among variables, there is an increasing demand for a more 
general theoretical framework that combines moderation and mediation analyses. Recently, 
statistical analysis of mediated moderation (meMO) effects has become a powerful tool for scientists 
to investigate complex processes. However, the traditional meMO model is formulated based on the 
homoscedasticity assumption, which is most likely to be violated when moderation effects exist. In 
addition, routinely reporting effect sizes has been recommended as the primary solution to the issue 
of overemphasis on significance testing. Appropriate effect sizes (ES) for measuring meMO effects 
are very important in reporting and interpreting inferential results. However, there does not exist an 
effective measure that allows us to answer the question regarding the extent to which a variable Z 
moderates the effect of X on Y via the mediator variable (M) in the meMO model. 

The article is organized as follows. First, the two-level moderated regression model proposed 
by Yuan, Cheng, & Maxwell (2014) was extended to a two-level mediated moderation (2meMO) 
model with single-level data, the statistical path diagram was structured according to the conceptual 
model and the equations. Second, several effect sizes were developed for the 2meMO effect by 
decomposing the total variance of the moderation effect. Third, to estimate the parameters of the 
2meMO model and the ES measures of the meMO effects, we developed a Bayesian estimation 


method to estimate the parameters of the 2meMO model. Fourth, a Monte Carlo simulation study 
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was conducted to evaluate the performance of the 2meMO model and the proposed ES measures 
against those with the meMO model. Finally, we illustrate the application of the new model and 
measures with a real data example. 

The simulation results indicate that the size of bias and MSE for parameter estimates are small 
under both meMO and 2meMO models whether the homoscedasticity assumption hold or not. The 
results of the coverage rate of the 95% CI for difmomg following 2meMO is comparable to those 
following meMO when the variance of moderation error is zero, which is the assumption the meMO 
model is based. However, when the moderation-error variance is nonzero, 2meMO yields more 
accurate estimates for difmemo than meMO does, the advantages of 2meMO over meMO become 
more obvious as the moderation-error variance increases. The results of Type I error rate indicate 
that 2meMO controls Type I error rather well, and the rates are close to 0.05 or below 0.05 under 
all the conditions. However, the Type I error rates of meMO tend to be higher than 0.05 when the 
moderation-error variance is nonzero. The power rates following the meMO and 2meMO models 
are comparable for the medium or large sample size, or when there is a large difference in meMO 
effects. While the value of power following 2meMO is slightly lower than that following meMO at 
small sample size, this result is mostly due to the inflated Type I error rate of meMO, and larger 
sample sizes and the smaller moderation-error variances correspond to more accurate estimates of 


P) 
mem 0: 


The results also indicate that, when the homoscedasticity assumption of the meMO model 
is satisfied, the effect size estimates following the two models are about the same. However, when 
the moderation-error variance is not zero, the results following 2meMO are more accurate than those 
following meMO. 

In summary, the article developed a 2meMO model with single-level data and proposed several 
measures to evaluate the size of the meMO effect explained by moderator variables in total, directly, 
or indirectly. The performance of the 2meMO model is compared against that of the traditional 
meMO model via Monte Carlo simulations. Results indicate that, when the assumption of 
homoscedasticity holds, 2meMO yields comparable results with those under meMO. When the 
homoscedasticity assumption is violated, estimates under 2meMO are more accurate than those 


under meMO. More importantly, the measures of the size of the meMO effect proposed in this article 


can be used as a supplement to the test of meMO effects and will meet the needs for reporting ES 
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in practice. Consequently, the 2meMO model is recommended for the analysis of mediated 
moderation, and the effect sizes (ESs) for the interpretation of the effect according to the questions 
of interest are better reported. 


Key words mediated moderation; heteroscedasticity; effect size; Bayesian estimation 
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